Quantum phase transition in a gapped Anderson model: NRG study 
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We use the numerical renormalization group method to investigate the spectral properties of 
a single-impurity Anderson model with a gap S across the Fermi level in the conduction-electron 
spectrum. For any finite 8 > 0, at half filling the ground state of the system is always a doublet. 
Away from half filling a quantum phase transition (QPT) occurs as function of the gap value S, 
and the system evolves from the strong-coupling (SC) Kondo - type state, corresponding to S < Sc 
toward a localized moment (LM) regime for 8 > Sc- The opening of the gap leads to the formation 
of one (two) bound states when the system is in the SC (LM) regime. The evolution across the 
QPT of their positions and the corresponding weights together with the dynamic properties of the 
model are investigated. 

PACS numbers: 71.15.Pd, 71.30.+h, 71.90.+q 



o 
o 



> 

(N 
(N 
O 

o 
o 



X 



I. INTRODUCTION 

One of the hallmarks of the Kondo effect^ is the raise 
of a narrow resonance, at the Fermi level, in the spectral- 
density function of a magnetic impurity embedded into a 
metallic host. The width of the resonance is proportional 
to the so-called Kondo temperature Tk, which is the 
characteristic energy scale. Below Tk the impurity spin 
is completely screened into a singlet by the host material. 
Below the Kondo scale (T < Tk), the low-temperature 
properties, such as the resistivity, spin susceptibility, or 
specific heat are properly described in terms of Landau 
theory 2 - of the Fermi liquid. 

The simplest approach to capture the Kondo effect is 
through the Anderson model4 It was used initially to de- 
scribe the formation of the localized moments (LMs) in 
metallic hosts. The model has inspired a lot of theoret- 
ical work, and a multitude of analytical and numerical 
methods was developed^. One of them, the numerical 
renormalization groupie— (NRG), originally proposed by 
Wilson, is known as one of the most reliable and accurate 
approach to capture the low temperature, low-energy 
physics of the model. Later, with the increase in the com- 
puting power, NRG was successfully extended to a broad 
range of more exotic quantum impurity models 7 - such 
as the two-channel Kondo problem^ coupled magnetic 
moments^ the coupling to a superconducting host^ or 
the soft-gap models 

In the present work we address a slightly different prob- 
lem. That of a magnetic impurity in a degenerate semi- 
conductor host which presents a gap across the Fermi 
level in the conduction-electron spectrum. In the nor- 
mal Anderson model the conduction band has a flat den- 
sity of states (DOS) at the Fermi level and the Kondo 
temperature is the only energy scale of the problem. 
The opening of a gap 8 in the conduction-band spec- 
trum introduces a new energy scale. The first ques- 
tion that arises is whether the Kondo state will survive? 
The problem was originally addressed by using different 



techniques: the quantum Monte Carlo^ density-matrix 
(DM) renormalization group^ Poor Man's scaling, and 
1 /N expansion^ with no consensus reached. Later, Chen 
and Jayaprakashi^ have used NRG method for the same 
problem. It was found that at half filling, any gap 6 > 
changes the ground state to a doublet. Away from half 
filling and for large enough gaps the Kondo state does not 
survive. More recently, the results obtained within the 
NRG framework were confirmed by using a local-moment 
approach.— 

The main goal of the present work is to extend the 
analysis of the model and to investigate its dynamic prop- 
erties and their evolution across the quantum phase tran- 
sition (QPT). In Sec. |TT]we present the theoretical model, 
then the results for the spectral properties are presented 
in Sec IIIII We give the conclusions in Sec IIVI 



II. MODEL AND NUMERICAL APPROACH 

To describe a local, quantum impurity state, coupled 
to a conduction band we use the generic Anderson model 
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Here £k is the host band dispersion, which is treated as 
a noninteracting one, td is the impurity-level energy, U is 
the on-site Coulomb energy at the impurity site, and V 
is the hybridization-matrix clement of the local impurity 
orbitals with the band states, which, in the present ap- 
proach, is considered momentum and spin independent. 
The number operator rid a = d\,d a , describes the occu- 
pation of the impurity level for spin-cr electrons. The 
mixing of impurity level with the host states is gener- 
ically described by the hybridization function A (u>) = 
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A R (uj) + i Ai (w) = l^ 2 Ek[ w - £ k+i V sgn (w)] . In 
general Ai(w) can be related to the density of states of 
the host band: Ai(w) = ttV 2 q(uj). In our model, a gap 
S > is present in the density of states. Then, Aj(ui) has 
the form 



A! (oj)=TQ(\u)\-6)e(D-\oj\) 



(2) 



with 2D the bandwidth of the host band. The normal 
Anderson model with a flat density of states is recov- 
ered in the limit of zero gap S — > 0, in which case the 
properly normalized density of states is q(ui) = 1/2D 
for uj G [— D, D] and the broadening function becomes 
r = nV 2 /2D. In the followings we will consider D as the 
energy unit. The real part of the hybridization function 
is obtained through the Hilbcrt transform as 
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The region of relevance corresponds to small gap val- 
ues, of the order of Kondo temperature ( we use the Hal- 
dane's expression for the Kondo temperature^ so it is 
properly defined only in the limit of zero gap, 6 = 0). 
We are mostly interested in the dynamical properties 
of the model. These are best described by the single- 
particle, retarded Green's function of the impurity site: 
G ret {t - t') = -i({d a (t),4(t')})e(t - t'). In gen- 
eral, the time-ordered Green's function of the interact- 
ing problem is given by the Dyson equation: G(oj) = 
[Gq(oj) — E(w)] -1 , in terms of the non-interacting, (U = 
0) Green's function Gq(uj) = [cj — td — A(w)] 1 and of 
the self-energy Then, the spectral representation 

A(uj) = —l/it^mGre^uj) can be readily obtained. 

Away from the half filling it can be readily shown,— by 
the simple perturbation theory at the Hartree-Fock level, 
that outside the gap, (|u| > 5) both real and imaginary 
parts of the self-energy are non-zero, so the spectrum is 
continuum. On the other hand inside the gap, (\uj\ < 5) 
the imaginary part of the self-energy vanishes. Because 
of that, the d- level Green's function has a single pole in- 
side the gap. This pole corresponds to a resonant state 
with a life-time inverse proportional with the imaginary 
part of the self-energy. Therefore, it corresponds to a real 
bound state with infinite lifetime. The position of the 
bound state (Ef,) can be obtained by solving the equa- 
tion G~e t (Eb) = for Ef,. At the non-interacting level, 
(neglecting the self-energy corrections) the bound-state 
energy is determined from 



Eb — ed = AR(£7b 
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Away from the half-filling this equation has always only 
one solution. When the self-energy corrections are in- 
cluded, it only leads to a renormalization of the local 
energy = + 5fte ~E(E{,), otherwise, the same scenario 




FIG. 1: Renormalization group flow diagram at T = for the 
gapped Anderson model, when the system is away from half 
filling. When 8 < 8c the SC fixed point is stable while for 
8 > 8c the LM becomes the stable fixed point. At half filling 
8c — so LM is stable irrespective of the value 8 > 0. 



holds, and still a single bound state is formed in the gap, 
but with the position slightly shifted. 

At half filling a more careful analysis is necessary. At 
the non-interacting level we do expect a single bound 
state to appear exactly at zero energy. However, the same 
perturbative analysis gives for the imaginary part of the 
self-energy inside the gap: 3mS(w) oc S(uj) and by the 
Hilbcrt transform the real part becomes 3?eS(o;) cx 1/oj. 
Using now Eq. Q, but with the self-energy correction 
included, it can be shown that at half filling, a pair of 
two, particle-hole symmetric poles are formed inside the 
gap. 

The perturbative approach described above provides 
us with helpful information relative to the formation of 
the bound states. Regardless of its simplicity it has 
its own limitations. Computing the self-energy in the 
second-order perturbation theory in U, the spectral func- 
tion of the impurity site can be resolved only at a qual- 
itative level. More than that, in the large-C/ limit the 
perturbation theory is supposed to fail. The previous 
analysis indicates that there is a fundamental difference 
between symmetric and asymmetric cases, but we can 
only speculate that the particle-symmetric case cannot 
be perturbatively connected to the non-interacting limit 
while the asymmetric model docs. At the same time no 
information relative to the nature of the ground state can 
be extracted. 

A much more rigorous analysis of the problem is pos- 
sible by using the NRG method^ which is a reliable 
approach to study a variety of quantum impurity mod- 
els. It consists of a logarithmic discretization of the host 
band into intervals of the form [A - ( n+1 ', A~ n ] with A 
some positive number larger than 1, (usually A ~ 2) and 
a mapping of the original Hamiltonian (fTJ) into a one- 
dimensional tight-binding chain (Wilson chain) such that 
the hopping couplings between nearest neighbors acquire 
a A dependence of the form cx A - ™/ 2 , followed by an it- 
erative diagonalization of the chain with one extra site 
added at each iteration step. In general the on-site ener- 
gies in the Wilson chain are also non-zero, but it can be 



analytically shown that for a host density of states with 
electron-hole symmetry they do vanish. The continuous 
limit corresponds to A — > 1 and for any A > 1 the NRG 
is an approximation. Keeping A small (~ 1.5) the com- 
puting time increases a lot while having a large A (~ 3.0) 
the accuracy of the calculation is compromised, especially 
for spectral properties at large energies. Therefore, in the 
present work, we present results for A = 2. There is ma- 
jor difference in the way the NRG is used when a gap is 
opened in the density of states. For example, at T = 0, 
in the normal Anderson model it is, in principle, possi- 
ble to increase the number of iterations to any value, but 
usually an upper limit is chosen such that, once the fixed 
point is reached, the NRG stops after a few iterations. 
In principle, the more iterations we use, the better the 
spectral functions at the Fermi level are resolved. On the 
other hand, for a gapped Anderson model the threshold 
5 fixes somehow the maximum number of iterations. Be- 
cause there are no longer states in the host band below 
5, we need to stop the NRG procedure at a given itera- 
tion N max which is gap dependent: N max = N{5) such 
that the typical energy scale A - ^""** -1 )' 2 is not much 
smaller than the gap 6. As a technical detail, we can take 
advantage of the symmetries of the model and classify the 
cigenstates of the Hamiltonian into multiplcts. Here, we 
have used two quantum numbers to label the multiplcts: 
(i) Q - the number of particle measured relative to the 
one particle per site; (ii) S - the total spin. In the normal 
Anderson model the relevant energy scale is the Kondo 
temperature. It can be defined through the Haldane's 
expression^ T% = y/UVjl e ^ d (e d +t/)/2t/r Thc gap s 

introduces a new energy scale. As we will see next, the 
relevant, dimensionless parameter that characterize the 
QPT is J/Tg.. 



III. SPECTRAL FUNCTIONS AND THE PHASE 
DIAGRAM 

We will start with a qualitative description of the 
renormalization group flow diagrams. These are dis- 
played in Fig. [TJ We will consider first the symmetrical 
model. For a normal Anderson model the strong-coupling 
(SC) fixed point is the only stable one, and the ground 
state is always a degenerate Kondo singlet characterized 
by the quantum numbers {Q, S) = (±1, 0). When a gap 
is opened in the density of states, the flow diagram com- 
pletely changes. The SC fixed point becomes unstable 
and the flow is toward the LM fixed point. At the same 
time the ground state changes to a doublet - (0, i), ir- 
respective of the value of the gap, so for the symmetric 
model the critical gap that describes the transition is 
6c = 0. 

Away from the half filling, and for the normal Ander- 
son model {5 = 0), probably the most interesting limit 
corresponds to the case when r <C — <C U. Corre- 
spondingly, the flow is toward the frozen impurity (FI) 
fixed point, which on the other hand can be identified 
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FIG. 2: (Color online) Spectral function for the symmetri- 
cal model for different gap values. The parameters used are 
-e<j = § = 0.2, T = 0.04. Tjj is the Kondo temperature in 
the absence of the gap. The up arrows indicate the positions 
of the bound states in each case, their magnitude being pro- 
portional to the corresponding weight rescaled with the gap 
value, i.e., TrTWb/S. Only the positive frequency is plotted, 
the spectrum for negative energies being symmetric. 
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FIG. 3: (Color online) The positions of the bound states and 
the corresponding spectral weights rescaled with the Kondo 
temperature T% (top layer) and with the gap value 5 (bottom 
layer) 



with the SC fixed point £ When a gap starts to open 
in the density of states the FI remains stable as long 
as the gap is smaller than a critical value 8c- For gap 
values larger than 5c {5c depends on the asymmetry of 
the problem, see Fig. 2]) the FI becomes unstable and 
the flow is toward the LM regime. At the same time 
the ground state changes accordingly as in the case of 
symmetrical model from a singlet to a doublet. 

In the following we will present results for the impurity 
spectral function, in these regimes. In general, within 
the NRG framework, the spectral function is given as a 
weighted sumi of 8 functions of the form 



A{uj) = Wrf (u - u>i) . 



(5) 
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FIG. 4: Phase diagram indicating the separation between SC 
and LM regimes. The parameter r\ — 1 + 2td/U describes the 
asymmetry of the system. At half filling, the critical gap is 
Sc — 0, and the system is always in the LM regime. The solid 
lines correspond to the analytical results obtained within the 
local-moment approach (see Ref. with similar rescaled 

coupling U — U/ttF. The inset displays the shift of the phase 
boundary line as function of A for 0.5 < rj < 0.75. 
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FIG. 5: Spectral functions for the asymmetrical case, close 
to the Fermi energy. In the upper panel the system is in 
the strong-coupling regime, characterized by 5 <g T%. In the 
lower panel the system is in the local-moment regime with the 
gap 8 of the order of T^- or larger. The arrows indicate the 
positions of the bound states and their amplitude is rescaled 
to nFWb/S. In the SC regime a single resonance develops 
while in the LM regime there are two. The dotted line in 
the left upper panel is the spectral function for the normal 
Anderson model. 



where the weights Wi can be computed directly with the 
NRG. To get a smooth spectral function the 8 functions 
need to be replaced by some smooth kernels such that the 
spectral sum rule J duiA(uj) = 1 remains valid. When in- 
vestigating the spectral properties of the gapped model 
we have used a slightly modified broadening procedure 
such that only the delta peaks which correspond to en- 
ergies outside the gap (\u>\ > 8) are broadened while for 
energies inside the gap (|w| < 8) the weights Wb and the 
corresponding bound-state energies Eb are extracted di- 
rectly from the excitation spectrum, so the spectral func- 
tion becomes 

A{u>) = A c ^(w)\\ u>s + Y^W h 8{u-E b )\ M<s . (6) 

b 

We will discuss separately the results for the cases 
when the system is at/away from half filling. In Fig. 
[2] we present typical results for the spectral function of 
the d-level of a symmetrical model. The presence of the 
gap in DOS preserves the electron- hole symmetry, so the 
spectrum remains symmetric. The black solid line is the 
spectral function for the normal Anderson model which 
develops the usual Kondo resonance below T^- Below 
the Kondo temperature the localized spin is screened by 
the conduction electrons, and the ground state is a sin- 
glet corresponding to S = 0. The opening of a gap in 
the density of states changes the physics dramatically. 
First of all, because no states are longer available at the 
Fermi level below the gap edge 8, the localized spin is 
no longer Kondo quenched into a singlet, and the ground 
state changes to a doublet (0, ^). 

At the same time, the spectral function A(uj) develops 
a gap and vanishes for energies |w| < 8 while two symmet- 



ric bound states at ±Eb with exactly the same weights 
develop inside the gap (Eb < 8). Although the ground 
state changes from a singlet to a doublet for any 8 
there are some reminiscence features of the Kondo peak 
and the behavior at energies > 8 resembles that of 
the normal Anderson model. For small enough gap val- 
ues the bound states are deep inside the gap. Increasing 
5, the Ei, moves toward the band edge. At the same 
time there is some transfer of spectral weight from the 
bound states to the continuum states. In the limit of 
8 3> T^- the bound states merge with the continuum, 
their spectral weights becoming vanishingly small and it 
cannot be resolved any more. In Fig. [2] the positions of 
the bound states are indicated by up arrows while their 
magnitudes are rescaled with the gap value: irTWb/S. 
In Fig. |3] the evolution of the rescaled bound-state en- 
ergy and weights as function of the gap is plotted in two 
different ways: (i) rescaled with Tj ( - the Kondo tem- 
perature; (ii) rescaled with 8 - the value of the gap. We 
have found that the energies Eb and the corresponding 
rescaled weights Wb = ir T Wb satisfy the following scaling 
equations in the limit S -C Tj^ : 

Eb_ f_S_V 



This scaling behavior can be understood in terms of 
an effective Hamiltonian. First the Anderson Hamilto- 
nian ([T]) is mapped by the help of the Schricffcr- Wolff 
transformations^ into a Kondo problem and then a low- 
energy, effective Hamiltonian (at the energy scale 8) is 
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FIG. 6: Typical evolution of the bound-state energies (left 
panel) and the corresponding weights (right panel) across the 
SC-LM quantum phase transition. In the left panel, the cir- 
cles (squares) represent the positions of the bound states in 
the SC(LM) regime. The crosses in the right panel indicate 
the sum of the weights of two bound states formed in the 
LM regime. The arrows point to the critical gap Sc where 
transition occurs. 



constructed, in which the impurity spin is coupled to the 
lowest electron/hole levels. The Hamiltonian is described 
in terms of an effective exchange coupling J e ff oc 
and a potential- like scattering term K e ff oc KqT^. In 
the presence of electron-hole symmetry, and in the LM 
regime, the scattering term vanishes: K e ff = 0, while 
only the exchange term survives. For a finite gap S the 
system is always in a LM state with a ground state in 
the (0, i) sector while the first excited state is in the 
(±1,0) sector. The energy Eh of the bound stated cor- 
responds to the transition between the lowest energies 
within these sectors [ (0, \) f-> (±1, 0) ] and this energy 



difference scales as Eh ~ Ei 
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Next, we will describe the quantum phase transition 
when the system is away from the half filling. Here the 
SC regime is much robust as compared to the symmetric 
case and the transition occurs always at a finite critical 
gap, 6c, such that for S < Sc the system is in the SC 
regime and for S > Sc the flow is toward the LM fixed 
point. The critical gap Sc depends on the asymmetry 
parameter 77 = 1 + 2ed/U. In Fig. @] we present the 
separation between these two regimes in the parameter 
space (77, 6c/Tft). We represent our numerical results 
for the critical values (symbols) together with the ones 
obtained by using the local-moment approach^ (solid 
lines). Close to the half filling the agreement between 
the methods is almost perfect, and only for 77 > 0.6 — 0.7 
deviations start to appear. The method of extracting 
the critical gap in the present approach is to some ex- 
tend different from the one used in the soft-gap Ander- 
son modelii where the phase transition boundaries can 
be extracted directly from the NRG flow diagram. This 
kind of analysis is not possible in our approach since the 
NRG gets truncated at the gap edge. In our approach the 
phase boundary was obtained by changing the asymme- 
try parameter for a fixed gap and counting the number 
of bound states in the gap. We believe that there is no 
other more efficient procedure for constructing the phase 
diagram for this model. More exactly, we have fixed the 



gap value S as well as the Coulomb interaction U, and we 
have changed the asymmetry parameter 77 = 1 + 2e c ;/C/. 
When the system is in the SC regime, away from the half 
filling, there is always one bound state in the gap, but 
decreasing the asymmetry, at some point we cross the 
phase boundary and the system evolves toward the LM 
fixed point with two bound states in the gap. To give a 
quantitative description of the transition, in terms of the 
spectral properties, when the gap is changed smoothly 
from 6 = to some large, S ^> T^, value, we have fo- 
cused on the regime with U = 0.4, = —0.05, and 
r = 0.01. For this set of parameters the Kondo tempera- 
ture is Tjf = 4.6 x 10~ 5 and the critical gap is found nu- 
merically to be Sc = 0.045 T^. In this particular regime 
(corresponding to T <C — <C U ) and for the normal 
Anderson model, below the Kondo temperature, the flow 
is always toward the SC fixed point, and the ground state 
is a non-degenerate singlet within the (—1,0) channel. If 
we steadily increase S, when S <C T^-, the ground state 
remains a singlet. At the same time a gap is opening in 
the spectral function A(uj) below energies \uj\ < 5. In the 
high-energy region, \uj\ > S, A(uj) has a similar structure 
with that corresponding to the normal Anderson model. 
The Kondo resonance forms below while the Hub- 
bard side peaks develop at similar energies. In this SC 
regime only one bound state develops inside the gap (top 
layer in Fig. [5]) in agreement with the perturbation the- 
ory. As approaching the transition S — > 5^ the energy 
of the bound state is shifted slowly toward the Fermi en- 
ergy (Eh —> 0) (see Fig. When <5 becomes larger than 
Sc the SC is no longer the fixed point, and the flow is 
toward the LM regime. The ground state changes to a 
doublet: (0, ^) while two non-symmetrical bound states 
develop in the gap. In Fig. |6]we present the evolution of 
the bound-states energy together with their characteris- 
tic weights across the transition point. On the LM side 
of the transition the localized states start to loose their 
weights as the gap is increased and their energy slowly 
merges into the continuum. 



IV. CONCLUSIONS 

Numerical renormalization group is by now a well- 
established method for studying correlation effects in 
quantum impurity models. We have applied it here to 
investigate the Anderson model with a gap in the conduc- 
tion band. By using a slightly modified version we were 
able to capture the modifications in the spectral prop- 
erties of the local operators across the quantum phase 
transition and the positions of the bound states formed 
inside the gap. In this way we were able to construct the 
phase diagram of the model. The NRG provides one of 
the most accurate tools for investigating the dynamical 
quantities of quantum impurity models. When DM- NRG 
is used (such as in our case) the sum rules for the spec- 
tral properties are satisfied up to the numerical precision. 
Therefore, the NRG results for the spectral functions are 
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net superior to those obtained within other approaches 
such as the LMA (see Ref. [l6| ) which do capture cor- 
rectly only the positions of the bound states while the 
spectral properties suffer because of the approximations 
made. 
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